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ABSTRACT 



The two and three-dimensional response of strong ocean currents to atmospheric forcing is 
studied using numerical simulations. In particular, surface cooling is explored as a possible 
mechanism for explaining an observed 100 km southward shift in the mean position of the Gulf 
Stream during winter. The magnitude and direction of the cross-stream circulation is highly 
dependent on whether or not a vertical mixing of momentum occurs when the water column 
convectively adjusts in response to the surface cooling. A weak cross-stream flow toward the 
higher sea-surface temperatures occurs in the surface layer if momentum mixing does not occur, 
whereas a stronger flow toward lower sea-surface temperatures results if momentum mixing does 
take place. The response of three-dimensional simulations is very similar to the two-dimensional 
simulations in the immediate vicinity of the front. The response due to horizontal cooling 
gradients is not large enough to displace the Gulf Stream appreciably southward in any of the 
numerical simulations. By contrast, a moderate increase in the zonal wind stress is more effective 
in displacing the core of the current system than are strong gradients in the surface cooling. The 
position of the Gulf Stream has also been observed to change as it encounters changes in the 
bottom topography. The simulated adjustment of flow to a 200 m seamount is consistent with 
the conservation of potential vorticity. The inclusion of surface forcing does not affect the 
adjustment to the topography in any of the simulations, and does not steer the flow toward (away 
from) the topography so that the characteristics of the downstream flow are changed. 
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I. INTRODUCTION 



Over the past decade there has been conflicting observational and theoretical evidence over the 
role of air-sea fluxes of sensible and latent heat during cyclogenesis. Using data from the Air Mass 
Transfer Experiment, Wei (1979) concluded that air- sea fluxes of sensible and latent heat did not 
contribute significantly to the available potential energy during cyclogenesis. However, there is 
evidence which suggests that ocean heat fluxes may be important for cyclogenesis at specific 
geographic locations. Baker (1979) finds that most cyclogenetic events in the southern hemisphere 
occur along the sub-polar ocean front. Sanders and Gyakum (1980) find that a positive 
correlation exists between areas of explosive cyclogenesis and locations of sea- surface temperature 
gradients. Sandgathe (1981) investigated changes in the evolution of model cyclones when the air- 
sea fluxes were removed. He selectively removed surface fluxes and found both changes in the 
intensification of an initial cyclone and in the development of secondary lows. The studies by 
Baker (1979), Sanders and Gyakum (1980) and Sandgathe (1981) suggest that knowledge of the 
locations of sea- surface temperature gradients may be needed for accurate prediction of 
cyclogenesis over the ocean. 

According to Sanders and Gyakum (1980), most explosive cyclogenetic events over the 
Atlantic Ocean occur in the vicinity of the Gulf Stream. There is a baroclinic zone associated with 
the Gulf Stream current which can be traced from Cape Hatteras eastward for hundreds of 
kilometers. The horizontal ocean temperature differences are not compensated by salinity effects, 
and there is an associated baroclinic jet with a front that behaves similarly to an atmospheric 
front. The local momentum balances in the Gulf Stream are similar to the momentum balances of 
an atmospheric front (Roden, 1976). The Gulf Stream temperature front tends to be vertical in 
the upper 50-100 m and then slopes steeply (about 100 m/10 km) toward the south with 
increasing depth. 
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Much study has been devoted to the understanding of the large-scale balances of the Gulf 
Stream. The review article by Veronis (1981) provides an excellent background on these studies. 
The investigation presented here deals with local balances and the effects of surface forcing on the 
Gulf Stream. The processes which form and maintain the Gulf Stream current system are not 
addressed. 

Most ocean front response studies have focused on wind forcing only (DeSzoeke, 1980; 
Cushman-Roisin, 1981; and Camerlengo, 1982). Recent studies have also considered the effects of 
a surface buoyancy flux in the forcing of oceanic fronts. Cushman-Roisin (1984) investigated the 
role of the large-scale buoyancy flux in maintaining the subtropical front in the North Pacific 
Ocean. Adamec and Garwood (1985) included a surface buoyancy flux in addition to wind forcing 
and demonstrated that the near-surface advective response of the ocean is enhanced when the 
surface buoyancy flux is downward (heating). 

Of particular interest is an investigation by Nof (1983) who studied the response of the Gulf 

Stream as it flowed into a region of cooling. The response of the Gulf Stream was examined 

because of the large cooling rates observed east of Cape Hatteras (Budyko, 1963; Worthington, 

1976). Nof assumed the density structure could be represented by two layers and that the lower 

-3 -3 

layer was motionless. A density increase of 1 x 10 gm cm over 1200 km m the downstream 
(eastward) direction was specified to simulate surface cooling of the ocean. Also, surface cooling 
was allowed to occur south of a 0.5 m s~ ^ eastward jet only. A southward cross-stream 
component developed in response to a zonal pressure gradient which resulted from the along- 
stream gradient in surface cooling (density increase). A power series solution for a steady state, 
inviscid, Boussinesq flow yielded cross-stream velocities of 3 cm s*^, which could lead to 
displacements as large as 90 km in a season. 

The displacements predicted by Nof appear to be consistent with Fuglister's (1972) isotherm 
analysis of the Gulf Stream during the 1965-1966 winter, which showed a 100 km north-to-south 
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displacement of the 5°C and 10°C isotherms. Halliwell and Mooers (1983) provide additional 
evidence for a southward shift in the Gulf Stream front during winter based on an empirical 
orthogonal function analysis of satellite infrared imagery. In that study, the average monthly 
position of long-period meanders is displaced 45 km southward between October and February, 
1975-1978. NoFs result that a steepening of the frontal interface will occur during winter is 
consistent with Worthington's (1976) isotherm analyses of the Gulf Stream east of the Grand 
Banks. 

The Gulf Stream is a highly variable meandering current system that parallels the temperature 
front. Because the front and the current penetrate very deep into the water column, the effect of 
gradients in the underlying topography may be important in determining the path of the Gulf 
Stream. The Gulf Stream frequently passes over areas where there are large changes in the bottom 
topography. The New England Seamount Chain in the western Atlantic has features which rise 3 
km above the surrounding ocean floor. There are observations which indicate that the seamount 
chain significantly affects the path of the Gulf Stream extension. Vastano and Warren (1976) 
observed drastic path changes as the Gulf Stream flowed over the New England Seamount Chain. 
Richardson (1981) observed that the New England Seamount Chain is a favored area for large- 
scale meanders and warm core ring generation. In addition, the topography appears to be 
responsible for a 100 km southeastward deflection of the current system. 

The overall objective of this research is to aid in understanding the effects of local 
atmospheric forcing on strong ocean currents such as the Gulf Stream. Particular attention is paid 
to those processes that displace the current system equatorward in the winter season. Because of a 
lack of detailed observations of the response of the Gulf Stream to surface forcing, and because 
the response of the ocean to surface forcing is expected to be nonlinear, numerical simulations are 
used to model the complex processes and response of the ocean. 
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The research is described in five major sections. Section II briefly describes the ocean models 
and initial conditions used in this investigation. The objective of Section III is to use numerical 
simulations to understand the effects of local atmospheric forcing on strong ocean currents. 
Processes that modify the oceanic thermal and momentum structure when a cross-stream cooling 
gradient is imposed are highlighted. A cross-stream cooling gradient is specified by extracting more 
heat on the warm side of the Gulf Stream front than on the cold side. The response due to an 
increase in the southward and eastward wind stress is also investigated to determine the relative 
importance of the additional effect of wind stress and those effects due to cooling only. Because 
downstream variations are not treated in these two-dimensional simulations, the Nof (1983) 
process is excluded. 

Section IV extends the results of Section III and includes the effects of along- stream variations 
in the forcing. The purpose of this section is to simulate numerically the three-dimensional 
response of an intense ocean current system such as the Gulf Stream as it flows into a region of 
surface cooling. The relative importance of the cross-stream and along-stream cooling gradients 
near the initial position of the associated baroclinic zone is discussed in detail. 

Section V uses a multi-level primitive equation model with a-coordinates to investigate the 
effect of wind forcing on a stratified flow similar to the Gulf Stream as it flows over bottom 
topography. One of the major goals of Section V is to determine if surface forcing affects the 
adjustment of the flow to the topography so that the mean position of the current system is 
significantly altered. Section VI summarizes the results of this research and suggests areas for 
future study. 
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II. THE MODEL 



A multi-level, nonlinear primitive equation model is used in these experiments to simulate the 
Gulf Stream response to atmospheric forcing. The governing equations in flux form for an ocean 
circulation model which predicts the shear currents in an ocean of uniform depth D are written in 



Cartesian coordinates as 

A 

U» = - ( uu )x - (vu) y - (wu) 2 + fv + A m V 4 u + K m u 22 + (uu) x + (vu) - (2.1) 

P 0 

/V 

V, = - (uv), - (w) - (wv) 2 - fu + A m V 4 v + K m v m + (uv) x + (vv) - ^ (2.2) 

P 0 

T t = - (uT) x - (vT) y - (wT) 2 + A t v*T + Kt-T^ (2.3) 

u* + v y + w 2 = 0 (2.4) 

P 2 = - Pg (2.5) 

P = Po (1 - a(T - T 0 )) (2.6) 



where ~ represents a vertical average over the entire depth and P is P - P, which can be calculated 
using (2.5). The symbols are defined in the List of Symbols. 

The ocean is assumed to be hydrostatic and incompressible, and density is prescribed to be a 
linear function of temperature alone. A rigid lid is specified by setting the vertical velocity to zero 
at the top and bottom of the domain, so that the integrated divergence over the entire depth of 
the water column is zero. Furthermore, the vertical average of the horizontal velocity components 
is zero, so that the model predicts only the vertical shear in the horizontal velocities. The 
barotropic components which are a result of a dynamic adjustment in the ocean are assumed to 
be smaller than the currents induced by the surface forcing. All calculations in the simulations 
presented throughout this study are performed on an f plane at 36°N. 
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The effects of diffusion are calculated using a biharmonic ( ) operator in the horizontal and 

2 2 

a Laplacian ( 3 / 3z ) operator in the vertical. The values for the horizontal eddy coefficients are 
chosen such that for the largest currents the Reynolds number based on the grid size is of order 1. 
That is, the advective and diffusive terms in the strong current regions will have a comparable 
magnitude for scales the order of Ax. The model also includes a generalized convective 
adjustment scheme based on a local gradient Richardson number: if the Richardson number falls 
below a critical value (0 for these simulations), then similar adjustments of heat and momentum 
occur so that the new adjusted profile will have a local gradient Richardson number slightly larger 
than the critical value. The details of the adjustment scheme are given in Adamec et al. (1981). 
Zero is used as the critical Richardson number in these simulations to isolate those effects due to 
convective overturning only in response to the surface cooling. The consequences of using a larger 
critical Richardson number are discussed in Section III. 

The response is due to surface forcing only, as there is no flux of heat or momentum across 
the north and south boundaries or the bottom. The surface boundary forcing due to wind stress 
and heat flux enters through boundary conditions on K ^ (u z , v z ) and Ky (T z ), respectively. 
Variations in the x-direction (downstream) are ignored for all two-dimensional simulations. The 
total horizontal extent (resolution) of the two-dimensional region is 400 km (2 km). The vertical 
resolution varies logarithmically from 6 m near the surface to 100 m at 1000 m in 20 levels. 

The horizontal dimensions (resolution) of the three-dimensional domain is 1280 km (20 km) 
in the east-west direction and 256 km (4 km) in the north-south direction. The vertical extent is 
1000 m, and the vertical resolution varies logarithmically from 10 m near the surface to 200 m at 
1000 m in 10 levels. The east and west lateral boundaries are open in the three-dimensional 
simulations and are prescribed using a variation of Orlanski's (1976) boundary condition as 
discussed by Camerlengo and O'Brien (1980). 
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The effect of topography is often included in baroclinic ocean models by specifying the vertical 
velocity as a bottom boundary condition, w = -V • VD, as in Huppert and Bryan (1975) and 
Semtner and Mintz (1977). The numerical model in these simulations uses the so-called a - 
coordinate transformation, which was introduced by Phillips (1957). The o -coordinate is a 
redefinition of the vertical coordinate that allows for variations in the depth of the fluid. 

The new vertical coordinate is defined by o = z/D, where D is the local depth of the fluid. 
The horizontal coordinate and velocity components are defined on constant o-surfaces. The 
governing equations for the flow are similar to the equations in Cartesian coordinates except for 
the pressure gradients. The horizontal gradient of pressure in o-coordinates is 

VP ~ % P ‘ V ° ’ ( 2-7) 

and the hydrostatic equation is 

P = -PgD. (2.8) 

The equations in a-coordinates are solved in a similar manner to that used for primitive equation 
models that use z as the vertical coordinate. The vertical resolution of sigma in these simulations 
varies logarithmically from 0.01 near the surface to 0.20 at a = -1 in 10 levels. The topography 
used is a Gaussian- shaped seamount which extends 200 m above a 1000 m ocean floor. The 
seamount extends approximately 120 km in the x and y directions. 

The initial conditions for temperature are derived from a smoothed version of a hydrographic 
section across the Gulf Stream at 38°N during mid-November 1970 (United States Naval 
Oceanographic Office, 1970) and are shown in Fig. 2-1. The horizontal sea-surface temperature 
gradient is 6° C in 20 km. On the northern side of the front, the temperature is uniform to a 

depth of 60 m, while the water is well-mixed to a depth of 100 m on the southern side. Notice 

also that the vertical stability below the mixed layer on the northern side is greater than the 
vertical stability on the southern side. As a result, the maximum horizontal temperature gradient 
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Figure 2-1. Initial temperature field of the upper 250 m in the 
simulated Gulf Stream front. The contour interval is 



occurs at a depth of 100 m. Assuming geostrophic balance, the initial zonal velocities (Fig. 2-2) 
are eastward with a maximum value at the surface near 2.0 m s" * . 

Rather than inferring a cooling rate from the observed along-stream density increase as in Nof 

(1983), a surface cooling is specified. The x-y dependence of the surface cooling (Fig. 2-2) is 

similar to the x-y dependence of the density increases used by Nof. The upward heat flux 

increases from zero at x = 320 km to a maximum value at x = 1000 km. Cooling is assumed to 

be directly proportional to the sea-surface temperature, so maximum cooling occurs on the warm 

side of the surface front. The initial cooling in the two-dimensional simulations is identical to the 

cooling in the three-dimensional simulations for x > 1000 km. As in Nof (1983), there is no 

-2 

cooling north of the front. The maximum value for the cooling is 970 W m" . 

The evaporative and sensible heat fluxes are directly proportional to the air-sea temperature 
difference in the bulk aerodynamic formulation. If the temperature of the air is assumed to be 
uniform over areas of non-uniform sea-surface temperatures, then the cooling will be largest 
where the sea-surface temperatures are highest. The choice of a cooling function which increases 
toward higher sea-surface temperatures seems justified. 

The choice of a cooling function which increases in the along-stream direction is more difficult 
to justify. The Air Mass Transformation Experiment (AMTEX) was conducted south of Japan 
during February 1974 and 1975. One of the scientific objectives of AMTEX was to investigate 
the physical processes by which cold continental air is modified through interactions with the sea. 
The experiment was conducted in a location where cold continental air originating from the Asian 
land mass was modified through interaction with the warm waters of the Kuroshio. An array of 
observing stations enabled a comprehensive study of the heat and momentum transfer between 
the atmosphere and the ocean. 

Kondo (1976) described the time evolution of the observed synoptic weather pattern, surface 
wind stress and sensible and latent heat fluxes during a cold outbreak (Fig. 2-3). As the cold air 
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Figure 2-2. Cross-section of the initial geostrophic velocities of 

the upper 250 m in the vicinity of the front (side face), 
and the east-west distribution of the surface cooling 
(top face). The contour interval for the velocities is 
40 cm s" , and the contour interval for the cooling is 
190 W m'2 (increasing toward the southeast). 
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Figure 2-3. Surface pressure field and the daily mean distribution 
of the surface wind stress (N-m), sensible heat flux 
(II m- 2 ) and latent heat flux (W m-2) during AMTFX 1974 
(from Kondo, 1976). The top row is February 22, the 
middle row is February 24 and the bottom row is 
February 26. 
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penetrates southward, the values of the surface heat and momentum fluxes become larger so that 
by 26 February the maximum wind stress is about 0.5 N-m (directed southward as seen in left 
frames) and the sensible and latent heat fluxes are in excess of 300 W m' and 800 W m , 
respectively. The maximum values of the heat and momentum fluxes occur initially near the 
coast, but progress further southward with the leading edge of the cold air. The maximum heat 
fluxes are located over the Kuroshio as the cold air is modified by interaction with the ocean 
(Kondo, 1976). An enlargement of the latent heat flux and surface wind stress patterns relative to 
the position of the Kuroshio on 26 February is shown in Fig. 2-4. The maximum values of the 
heat and momentum fluxes occur to the south of the mean path of the Kuroshio, and there is a 
cross-stream gradient in the heat flux. Notice that the Kuroshio is approaching these maxima 
from the south, so that the cooling and wind stress is increasing in the downstream direction. 
The specification of cooling increases in the downstream direction in the numerical simulations is 
consistent with these measurements. 

Observations in the Atlantic support the choice of a cooling function which increases in the 

downstream direction for seasonal time scales. Gorshkov (1978) presents a monthly average of 

latent heat flux over the North Atlantic during February (Fig. 2-5) based on 70 years of ship 

observations. The maximum values occur off the east coast of the United States and are 

concurrent with the mean position of the Gulf Stream. The Gulf Stream flows from the south 

toward these maximum values and thus experiences a downstream increase in surface cooling. 

-2 

The maximum value of the average latent heat flux is 180 W m~ . 

In all of the simulations presented in this investigation, the surface forcing is applied 

impulsively. A time series of the computed latent heat flux from a study by Agee and Howley 

(1977) at Station Ryofu, which was located within the path of the Kuroshio during AMTEX 

1974, is shown in Fig. 2-6. A rapid increase in surface cooling is evident at the initial stages of the 

-2 

cold outbreak on 22 February. The values of the latent heat flux rise from 50 W m to over 800 
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February 

o 

Figure 2-6. Time series of latent heat flux (W m” ) at Station Ryofu during 
AMTEX 1974 (from Agee and Howley, 1977). 
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W m . The computed value of the latent heat flux combined with an estimated maximum 

-2 

sensible heat flux of 300 W m is about 15% greater than the maximum heat flux used in the 
simulations (970 W m ). Similar cooling rates during cold outbreaks have also been observed in 
the Gulf of Mexico by Huh et al. (1984). 

These simulations also investigate the role of wind forcing as a possible mechanism for shifting 
the Gulf Stream. The maximum wind stress used in these simulations is 0.2 N-m, w r hich is near 
the low end of values for the observed wind stress during the cold outbreak during AMTEX 1974 
(Fig. 2-4). Two directions of the wind stress are considered. The first is from the north, w T hich is 
the dominant direction of the wind during the outbreak of 22 February 1974 (Fig. 2-3), and the 
other is from the west. An increase in the eastward wind stress is considered because Gorshkov 
(1978) found that the average wind over the Gulf Stream increases to 5-7 m s’* in the eastward 
direction during winter. A value of 0.2 N-m is consistent with the observed values of the wind 
stress. The simulations presented here are designed to study the response of the Gulf Stream once 
it has left the coast near Cape Hatteras (Fig. 2-5). The observations during AMTEX 1974 were 
taken over the northward flowing Kuroshio before it leaves the coast of Asia. It is felt that the 
relative direction of the applied wind stress to the Gulf Stream realistically depicts outbreaks 
which originate from the cold continent to the north and west. 
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III. TWO-DIMENSIONAL SIMULATIONS 



The experiments in this section examine the two-dimensional effect of a horizontal gradient in 
the cooling across a strong ocean frontal feature such as the Gulf Stream. Such a gradient in 
cooling may occur when a nearly uniform mass of cold air flows over the sea-surface temperature 
gradient. An analytic model with very simple mixed-layer dynamics is used first to aid in 
understanding the response of an ocean front to cross-stream cooling. A more sophisticated 
numerical solution of the problem follows in later sections. 

Consider a hydrostatic fluid whose temperature is a constant T j down to depth z = -H, and 
then a constant T 2 for the remainder of the water column to depth z = -D. The governing 
equations for linear in viscid shear flow on an f plane are given by 

u^ — fv = 0 (3.1) 

/V 

v t + fu = - L (3.2) 

P 0 



where u and v are the Eulerian horizontal components of velocity, and p Q is the density of the 
fluid at a reference temperature T q . In this study, x is positive to the east and is the downstream 
direction, y is positive to the north and is the cross-stream direction and z is positive up. For 
simplicity, variations in the x-direction have been ignored and the wind stress is assumed to be 
zero. A single equation for v can be found by taking a local time derivative of (3.2) and 
substituting from (3.1) to obtain 

/V 

v„ + f J v = - (3.3) 

P o 



If density is prescribed to be a linear function of temperature alone, then it is possible to write an 

expression for P in terms of temperature as 

0 00 

P = J P g(l - a( T-TJ) dA - ± JJ Pg(l - a(T-T 0 )) dA dz (3.4) 



-D z 
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where a is the coefficient of thermal expansion and g is the acceleration due to gravity. The first 
term on the right side of (3.4) is the pressure, and the second term is the vertical average of the 
pressure. Taking a local time derivative of (3.4) removes all terms in the integrand which are not 
multiplied by T^. If changes in temperature are due solely to surface heat flux, which is assumed 
to penetrate to depth z = -H, then the local rate of change of temperature can be written 



T t = - 



T t = 0 



P„C p H 



0 > z > -H 



-H > z > -D 



(3.5) 



where Q is the upward surface heat flux. There is no change in temperature below z = -H. 
Substituting (3.5) into the time derivative of the first term on the right side of (3.4) yields 

0 > z > -H 



P. = - -2§- Q — 
‘ C p V H 



P t = -Pr- Q 



-H> z > -U 



(3.6) 



or by performing the integrations and removing the vertical average, 

0 > z > -H 



p = _ O (— + 1 - — ) 
‘ C p ^ H 2D J 

P t =^i Q JL 

1 C„ 2D 



-H > z > -D 



(3.7) 



Differentiating (3.7) with respect to y and substituting the expression for into (3.3) yields an 
equation for v forced by a gradient in the surface cooling 

0 > z > -H 



'• + fV =^ Q ’<F + '-£> 



V + f 2 v = 5S- Q — 

tt ^ v y or 



P c 

p 



2D 



-H > z > -D 



(3.8) 



It was assumed in (3.8) that H is not a function of y, which is not the case for the initial 
temperature structure in the numerical simulations (Fig. 2-1). A horizontal temperature gradient 
can also be induced by variations in H with uniform surface cooling. The contribution of the 
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variation in the upper layer thickness ( 5H) to a variation in the surface cooling ( <5Q) for forming 
horizontal temperature gradients is measured by the ratio ( <5H/H): ( SQ/Q). For the initial 
conditions and cooling used in the numerical simulations, ( 8 H/H) ~ 0.1, while ( 8 Q/Q) ~ 1.0. 
Thus, the ratio is about 0.1, which indicates that the gradient in the surface cooling is more 
important in forming horizontal temperature differences than is the gradient in the thickness of 
the upper layer. This justifies the derivation of (3.8) which will be used to interpret the numerical 
solutions shown in this section. 

The free response of (3.8) is inertial motion. Obtaining a particular solution to the forced 
problem requires assumptions about the heating rate Q. If the right side of (3.8) is assumed to be 
constant in time, then a particular steady state solution of (3.8) is 

v = - ag -- Q. C— + 1 — — ) 0 > z > -H 

P 0 C p 2D J - ~ 

Q, 1L -H > z > -D (3.9) 

P oW 

The forced solution in the upper layer is a cross-stream velocity that is directly proportional to 
and in the direction of the horizontal gradient of upward heat flux. 

An interesting application of this simple model can be made in the case of an existing surface 
temperature front. If a bulk aerodynamic formulation is assumed, then Q (latent and sensible heat 
fluxes) is directly proportional to the air-sea temperature difference. If the air temperature is 
constant, then Q has the same form as the sea-surface temperature. Consider now the special case 
of sea-surface temperature decreasing linearly with y, i.e. colder to the north, then < 0 is a 
constant. The forced solution near the surface is a southward cross-stream velocity. The physical 
interpretation is that a gradient in the surface cooling reduces the horizontal temperature gradient 
in the upper layer and consequently reduces the horizontal pressure gradient. If the front is 
initially in geostrophic balance, the cooling creates a supergeo strophic flow, i.e. the magnitude of 
the Coriolis term exceeds the presure gradient term. The result is a cross-stream current directed 
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to the right of the geostrophic flow, or in the direction of the larger upward surface heat flux. It is 

interesting to notice that the cross-stream velocity predicted by (3.9) is independent of time if the 

-2 3 1 

heat flux remains constant. Using values of 9.81 m s for g, 1000 kg m for p Q , 4186 J kg 

for Cp , 2 x 10 °C‘* for a, 100 m for H, 1000 m for D, and 970 W nT^ in 20 km for Q at 

36 N yields a value of 0.31 cm s~* for the surface cross-stream velocity. The average of the cross- 

stream velocity for the upper layer is 0.16 cm s'* 

It is possible to show that the linear forced response in the along-stream flow is given by 

u(t) = u(0) + A Q, (-jj + 1 - Jr) t o > z > -H 

PoCpf H 2D 

u(t) = u(0) - Q f ^-t -H > z > -D (3.10) 

which is linearly dependent on time. As long as a north-south gradient in the upward heat flux 
persists, the upper layer horizontal temperature and pressure gradients and associated zonal 
geostrophic velocities are diminished. 

The implication is that cross-stream cooling due to an atmospheric heat flux will both reduce 
and displace southward the north-south temperature gradient at the sea surface. While the 
induced cross-stream velocity is not very large (about 1/10 the value in the process treated by Nof, 
1983), additional effects not considered by the linear theory above, such as the vertical exchange 
of momentum due to convection, may have a significant effect on the internal response of the 
ocean. 

The numerical experiments designed to extend the above linear analysis are divided into three 
classes. The first class includes cases in which vertical redistribution of heat only is permitted as 
the upper layer cools and sinks due to interaction with the cold atmosphere. The second class 
includes both the exchange of heat and of momentum in the convective adjustment problem. 
That is, the convective overturning in the column is assumed to mix momentum as well as heat 
uniformly in the vertical. The third class includes cases which consider the additional effect of an 
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increase in the southward and eastward wind stress. Both heat and momentum are allowed to be 



mixed by the convective overturning when wind stress forcing is included. 

In each of the simulations in this section, the response fields will be shown as a deviation from 
a control simulation which uses zero forcing. Because of the large gradients in temperature and 
zonal velocity, part of the response in any simulation will be due to diffusive effects only. Since 
those effects are not the primary interest, the diffusive response in a zero -forcing simulation is 
removed to isolate the response due to surface forcing. In addition, the flow fields are averaged 
over 21 hours to filter the inertial response. Hence a flow field at hour 72 is actually the average of 
the flow for hours 51-72. 

A. EFFECT OF VERTICAL MIXING OF HEAT ONLY (EXPERIMENT 3-1) 

A southward (negative) surface v component is evident after 72 hours of integration as shown 
in Fig. 3-1. The response in the v field is mainly confined to the upper 100 m in the region of the 
surface front. The maximum magnitude is slightly greater than 0.3 cm s"^’ which is in good 
agreement with the above linear theory. The instantaneous values of the v components at hour 72 
(not shown) are remarkably similar to the averaged v components shown in Fig. 3-1, which 
indicates almost no time dependence in the forced flow. 

To illustrate the physical processes that determine the response of the front to the imposed 
cooling, the contribution of the individual terms to the time tendencies in (2.1), (2.2) and (2.3) 
were calculated. An inspection of the terms comprising the average v tendency for hours 51-72 
reveals that the difference between the pressure gradient and Coriolis terms is an order of 
magnitude larger than any of the other terms in (2.2). Thus, the tendency for the cross-stream 
flow is being accurately described by (3.2). As in the analytic work above, the surface cooling is 
reducing the north-south pressure gradient so that the Coriolis term is slightly larger than the 
pressure gradient term. The imbalance produces a southward flow opposite to the mean pressure 
gradient near the surface front. The along-front flow is constantly adjusted so that by hour 72, 
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Distance (km) 

Figure 3-1. Average of the v component of velocity for hours 51-72 for Experiment 3-1. 

The contour interval is 0.1 cm s - ^ and negative contours are dashed. The 
flow is positive northward (to the left). The arrow indicates the position 
of the warm edge of the surface front. 



the magnitude of the u component of velocity has been reduced by over 0.08 ms**’as shown in 
Fig. 3-2. The change in u is also confined to the upper 100 m in the region near the surface front. 

Although the sense of the horizontal cooling gradient is to destroy the ocean front, the 
dynamic response of the near-surface flow is thermally indirect (warm water sinking and cold 
water rising) and this tends to steepen the front and maintain it at its initial position. Notice that 
the largest southward flow occurs on the warm edge of the front where the cooling is greatest. 
Hence, the flow is divergent on the cool side of the front and upweiling occurs. There is also a 
small region of positive v near the surface at y = -80 km, which indicates convergence and 
downwelling on the warm side of the front. An indication of the effect of the indirect circulation is 
shown in Fig. 3-3 which illustrates the temperature tendency due to the combined effect of 
horizontal and vertical advection. Advective processes increase the temperature on the warm side 
of the front and lower the temperature on the cold side of the front. The zero tendency line is 
located at the warm edge of the front, and the strongest tendencies occur at the surface on either 
side of the warm edge of the front. The horizontal advective term is the contribution from the 
effect of confluence (difluence) in the strengthening (weakening) of the horizontal temperature 
gradient. The vertical advective term represents the effect of a differential tilting of the isotherms 
from the horizontal. 

By taking a horizontal derivative of the temperature tendency field, one obtains the 
frontogenetic tendency. The contribution from the advective terms is closely related to the 
geostrophic frontogenetic tendencies discussed by both Sawyer (1956) and Eliassen (1962). If 
frontogenesis is defined as areas where d/dt (- dT/ dy) > 0, then the advective tendency for 
fronto genesis is positive at the warm edge of the surface front and negative on either side of the 
warm edge. The nonlinear terms are acting to strengthen the surface front at the warm edge of the 
front. Sawyer and Eliassen also considered the effect of a cross-front shear on the horizontal 
temperature gradient. Since no variation in the x-direction is included here, the effect of cross- 
front shear is zero. 
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Figure 3-2. Average of the change in the u component of velocity for hours 51-72 for 
Experiment 3-1. The contour interval is 1 cm s"^ and negative contours 
are dashed. 
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Distance (1cm) 

Figure 3-3. Combined effect of horizontal and vertical advection on the temperature 

tendency at hour 72 for Experiment 3-1 The values are multiplied by 10. 
The contour interva' is 0.5 x 10"^ °C~’ and negative contours are dashed. 



The temperature change after 72 hours (Fig. 3-4) reflects an advective effect. A maximum 
temperature change of -0.82°C occurs south of the warm edge of the surface front. The 
temperature response closer to the warm edge of the surface front, but where the cooling is just as 
strong, is not as great (-0.69°C) due to the tendency of the horizontal convergence of heat to 
oppose the effect of surface cooling. For this simulation, the vertical extent of the cooling is 
constrained by the initial well-mixed layer. 

In summary, a horizontal cooling gradient produces a very small cross-stream velocity which is 
not likely to affect the position of a baroclinic zone as strong as the Gulf Stream. Although the 
cooling produces a southward component in the flow, any tendency to displace the surface front 
is offset by the vertical velocity tending to steepen the front and maintain the front at its initial 
position. The temperature response includes changes that are not uniform over areas where the 
surface heat flux is uniform. The largest temperature response occurs at the warm edge of the 
surface front. 

B. EFFECT OF MOMENTUM MIXING (EXPERIMENT 3-2) 

The convective adjustment used above exchanges heat between adjacent levels to produce a 
temperature profile that is at least neutrally buoyant. Unlike temperature, exchanges of 
momentum usually occur on the same time scales of the mean flow or in this case, 1/f. However, 
if surface cooling is forcing convective overturning, very little can be said about the time scale of 
the turbulent exchange of momentum. If the convection is sufficiently strong, then all properties 
of the water are likely to be mixed between adjacent layers. In laboratory experiments of flows 
with large Reynolds number, Thorpe (1973) found that both the heat and momentum profiles 
attained similar shapes when the Richardson number dropped below a critical value for short 
(order of seconds) time scales, which indicates that heat and momentum were being exchanged 
similarly. Although the exchange of momentum between layers has not been investigated in great 
detail, there have been numerical simulations which have shown the exchange of momentum in 
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0.25 °C and negative contours are dashed. 



an adjustment problem to be an important influence (Adamec el ai , 1981; Schopf and Cane, 
1983). 

The response in the frontal simulation with identical forcing to Experiment 3-1, but with 
momentum mixing included, is very different from the previous simulation. The v components 
after 72 hours for this simulation are shown in Fig. 3-5. The most important differences from 
Experiment 3-1 are the reversal of sign near the surface and that the values are an order of 
magnitude larger. Below the surface, the flow is to the south. Here, as in Experiment 3-1, the 
horizontal pressure gradient and Coriolis acceleration are the largest terms in the v momentum 
equation by at least one order of magnitude. However, the absolute value of the pressure gradient 
force is larger than the absolute value of the Coriolis acceleration near the surface front. 

The difference in the response of the v components can simply be explained in terms of the 
momentum mixing. When surface cooling occurs, the gradient Richardson number is less than the 
critical value. The relatively high momentum near the surface is mixed with relatively low 
momentum immediately below, so that the surface currents decrease relative to their initial values. 
Immediately below the surface layer, the currents increase. As the Coriolis term decreases at the 
surface, and the pressure gradient term is not decreased by a similar amount, the flow becomes 
sub-geostrophic and a flow down the mean pressure gradient occurs (positive v components, Fig. 
3-5). Below the surface layer where u has increased due to the turbulent exchange with relatively 
higher momentum from above, the Coriolis term is larger than the pressure gradient, and the 
response in v is up the mean pressure gradient (negative). 

The changes in the u component after 72 hours (Fig. 3-6) are consistent with these arguments. 
Notice that the speed change in the surface layer is negative (reduction of the jet) and twice as 
large as the changes in Experiment 3- 1 due to the additional effect of mixing with relatively lower 
momentum from below. Immediately below the surface layer, the current speed has increased 
relative to the initial value. At 60 m depth, there is an area with reduced speeds. This is a layer in 
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Figure 3-5. As in Fig. 3-1 except for Experiment 3-2 and the contour interval is 
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Figure 3-6. As in Fig. 3-2 except for Experiment 3-2 and the contour interval is 



which the initial vertical shear in u is a maximum. Above this area, the shear is smaller and 
constant with depth. Locations where the vertical shear is large are locations where momentum 
mixing will also be very large. Plots of the instantaneous currents reveal a high variability in the 
currents at depth. The variability appears to be inertial. There is little evidence for propagation of 
the inertial energy vertically or horizontally. The oscillation remains confined to the upper 100 m 
toward the warm edge of the front. 

The reversal in sign of the surface v components induces vertical velocities which tend to 
reduce the slope of the front. The vertical components of velocity at hour 72 are shown in Fig. 
3-7. Notice the upwelling on the warm side of the front where downwelling occurred in the 
previous case. On the cool side of the front, downwelling is taking place. The sense of the vertical 
advection is direct (warm water rising and cold water sinking). The divergence on the warm side 
of the front and the convergence on the cold side of the front tend to diminish the temperature 
gradient across the front. 

In the atmosphere, there have been several investigations that show that upper level fronts in 
the vicinity of jets are preferred regions for vertical mixing (e.g. Kennedy and Shapiro, 1980). 
Although large vertical eddy diffusivities are apparent in the upper atmosphere, care must be 
taken in applying those results to ocean surface fronts since upper-atmospheric fronts are far 
removed from the effects of a planetary boundary layer. In a numerical simulation of upper- 
atmospheric frontogenesis, Shapiro (1982) notes that if a maximum in the heat flux occurs in the 
vicinity of a front, the forced circulation is thermally indirect. However, if the turbulent 
momentum fluxes dominate the heat fluxes, the forced circulation will be thermally direct, which 
is consistent with the simulations in this study. Other studies have indicated a dominance of 
turbulent momentum exchanges over turbulent heat exchange. Gidel and Shapiro (1979) 
investigated the role of turbulent heat and momentum fluxes in the production of potential 
vorticity near an upper tropospheric front for a range of Prandtl numbers (K^/K^). For a 
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dashed. The values have been multiplied by 100. 



Prandtl number of 1, the change in vorticity due to the turbulent momentum flux dominated the 
change in vorticity due to the turbulent heat flux. 

Each of the above simulations was repeated with different cooling rates. The sense of the 
ocean response is not dependent upon the magnitude of the cooling, however the magnitude of 
the response is nearly linearly dependent on the magnitude of the cooling. When the cooling rate 
is reduced by 50%, the changes in the temperature and current fields are similarly reduced by 
about 50%. 

The simulations which included momentum mixing were repeated using a critical Richardson 
number of 0.25 to include possible dynamic instability effects. The response of the cross-stream 
and the along- stream velocity components are in the same direction as in the case with critical 
Richardson number equal to zero, but the magnitude of the currents is increased by about 30%. 
The increase in the response is due to relatively high momentum at the surface being mixed with 
relatively low momentum from below as the condition for instability (mixing) is satisfied more 
frequently. However, the large stratification in the thermocline restricts the effects due to 
momentum mixing to the upper layers. Although, the response in the currents is enhanced, the 
responses still occur in the same locations. 

The role of momentum mixing in the ocean is not well documented. Numerical simulations 
by Adamec et aL (1981) show that including mixing of momentum alters significantly the oceanic 
response due to hurricane forcing. Schopf and Cane (1983) found that entrainment of momentum 
is a non-negligible process in maintaining equatorial undercurrents. The role of momentum 
mixing is particularly interesting in the simulations of the response of the Gulf Stream to surface 
cooling since the mixing of momentum produces a very different response in the along-front flow 
than in the simulations without it. Measurements of the change in the Gulf Stream velocity 
profile during a cold outbreak are required to document the importance of this process in the 
ocean. 
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C. EFFECT OF WIND FORCING (EXPERIMENTS 3-3 AND 3-4) 



Two components of the total heat flux, sensible and latent heating, are directly proportional to 
wind speed in the bulk aerodynamic formulations. Large losses of heat by the ocean are normally 
accompanied by an increase in the wind speed. In the previous simulations that included the 
effects of cooling only, the treatment of the vertical structure of momentum was a critical factor in 
determining the magnitude and direction of the internal ocean response. In this section, a surface 
momentum flux is also included to determine if it is an important factor in determining the 
response of the ocean to a cold air outbreak. 

Two wind forcings are considered: a cross-stream southward wind stress simulating an 
outbreak of cold air from polar regions, and such as the cold outbreak during AMTEX 74 (Agee 
and Howley, 1976), and a wind stress in the direction of the along-stream flow to simulate a cold 
air outbreak from the cold continental area to the west. For each of these simulations, the 
maximum wind stress is 0.2 N-m and is applied impulsively. This wind stress represents the low 
end of values observed by Agee and Howley (1976). The effects of the combined cooling and 
wind forcing are compared to the effects due to cooling only. 

The change in u velocity for the simulation with cooling and a southward wind stress is shown 
in Fig. 3-8. The difference from the cooling-only simulation is a further decrease in the speed of 
the surface jet by 0.1 m s*' The reduction of the eastward component is consistent with an 
Ekman-type response to the southward wind stress. The reduction of the jet is evident 
throughout the surface layers as the convective overturning efficiently mixes momentum. Iffiere 
are areas below the warm edge of the surface front where the jet speeds have increased. Notice 
that the reduction in the jet in the cooling-only simulation is larger than the additional reduction 
due to the wind forcing alone (0.16 m s as opposed to the additional 0.1 m s ^). However, 
significant reductions in the initial jet penetrate to depths twice as large when the effect of wind 
forcing is included when compared to the cooling-only simulation. 
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Figure 3-8. As in Fig. 3-2 except for Experiment 3-3 and the contour interval is 



The effect of the southward wind stress on the v velocity is shown in Fig. 3-9. The large 
positive v components near the surface in the cooling-only simulation are greatly reduced. The 
tendency to produce a current toward the north due to momentum mixing (Fig. 3-5) is opposed 
by the effect of the wind stress. Farther to the north and south of the front, the v components are 
negative, which is in the direction of the applied wind stress. The response immediately below the 
surface front is also somewhat diminished due to the opposing tendencies of momentum mixing 
to produce a sub-surface current toward the south and the wind stress which forces a sub-surface 
current toward the north. 

On the warm edge of the front the flow is divergent, and toward the cooler side, the flow is 
convergent. The response in this simulation, as in the cooling-only simulation, is a thermally 
direct circulation at the warm edge of the front. The strength and shape of the vertical circulation 
pattern (Fig. 3-10) is similar to the cooling-only simulation. 

The response in the temperature field (not shown) is very' similar to the cooling-only 
simulation. There is a slight southward shift in the temperature response due to the advection 
induced by the southward wind stress. The temperature response, as in the previous simulations, 
is dominated by the surface heat flux and not by the internal response of the ocean. 

The response to the inclusion of an eastward wind stress is very dramatic. The sense of the 
Ekman advection for this simulation is southward. The changes in the near-surface v velocity at 
hour 72 (Fig. 3-11) are negative (southward) everywhere except for a small area at the surface. 
The maximum change occurs in the vicinity of the front where the v components are 0.04 m s" ^ 
In the cooling-only simulation, the maximum v components are in the same area and of the same 
magnitude, but in the opposite direction. Because the magnitude of the Ekman response is 
inversely proportional to mixed-layer depth, the surface v velocities are divergent at the front. In a 
similar simulation of the response of the Maltese ocean front to wind forcing, Adamec and 
Garwood (1985) found that the forced divergence was very effective in spreading the surface 
isotherms and reducing the intensity of the flow near the front. 
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Figure 3-9. As in Fig. 3-1 except for Experiment 3-3 and the contour interval is 
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Figure 3-10. As in Fig. 3-/ except for Experiment 3-3. 




I I I 



(m) in ds a 



50 



0 -40 -80 -120 -160 -200 

Distance (1cm) 

Figure 3-1 *1 . As in Fig. 3-1 except for Experiment 3-4, the contour interval is 1 cm s - ^ , 
and for values of y between -200 and 0 km. 



The u components for the simulation that included an eastward wind stress are shown in Fig. 
3-12. Near the initial position of the front, the u components are decreased by over 0.38 m s~^ 
whereas the u components are increased by 0.12 m s~^ farther to the south. This southward shift 
of the jet is due to the geostrophic adjustment of the flow to a southw T ard shift in the horizontal 
temperature gradient. Colder water from the north is advected southward and sinks as it 
encounters the warmer water in the south. The tendency for advection is to diminish the 
northern edge of the horizontal temperature gradient and form a relatively stronger temperature 
gradient to the south. As a result, the u components decrease on the northern side of the initial 
jet and increase to the south. In the region near -140 km, a southward advection causes positive 
(negative) changes in u where the initial du/ dy is negative (positive). 

The maximum surface temperature change (Fig. 3-13) at the front is more than 1.0°C lower 
than the temperature changes in the cooling-only simulation due to the southward advection of 
cold water. At y = -30 km and y = -140 km, there are areas of temperature increases below the 
initial isothermal layer. This warming is due to the effect of downwelling occurring in those areas. 
The enhanced cooling at y = -120 km is similarly due to a local area of upwelling. As in all 
cases that included momentum mixing, the vertical circulation near the front (not shown) is 
thermally direct. 

In summary, a moderate increase in the eastward wind stress has a greater effect on the 
position of the current system than does a strong cross-stream cooling gradient. Only the 
simulation which considers the response to an eastward wind stress has a noticeable shift in the 
position of the current system after three days of forcing. Nof (1983) considered a different 
mechanism to produce cross-stream currents. Nof assumed that an adjustment by the ocean to an 
along-stream cooling gradient produces a cross-stream current of sufficient magnitude to displace 
the Gulf Stream by as much as 90 km in a season. Nof s mechanism must be modelled in three 
dimensions. Since the direction and magnitude of the cross-stream current in the two-dimensional 
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Figure 3-12. As in Fig. 3-2 except for Experiment 3-4 and for values of y between -200 
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Figure 3-13. As in Fig. 3-4 except for Experiment 3-4 and for values of y between -200 



simulations are dependent on whether or not momentum mixing is included as part of the 
convective adjustment to surface cooling, it is expected that momentum mixing will also play an 
important role in three-dimensional simulations. 
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IV. THREE-DIMENSIONAL SIMULATIONS 



Three sets of experiments are designed to investigate the three-dimensional response of a 
strong ocean current to local atmospheric forcing. The purpose of the first two sets of experiments 
is to compare the relative importance along-stream cooling gradients to the cross-stream cooling 
gradients treated in Section III in the two-dimensional simulations. The design of these 
experiments is similar to the design of the experiments in Section III. In the first experiment, 
only the vertical mixing of heat in the convective adjustment problem is considered. The second 
experiment is similar to the first but also includes the effect of vertical mixing of momentum as 
part of convective adjustment process. The third set of experiments considers the response of the 
ocean currents to an increase in the zonal wind stress. Part of the response in the horizontal 
velocity components when momentum mixing is included is due to inertial motion. Although 
inertial motion is not of primary interest to this investigation, an interesting effect occurs which is 
due to advection. The strong zonal advection alters the period of the oscillations in these 
simulations so that the period is near 15 hours rather than the expected 21-hour inertial period. 
An analysis of how advection affects the period of oscillation in these simulations is presented in 
Appendix A. For all the results presented, the flow fields have been low-passed filtered to remove 
the inertial response, and the fields are presented as a deviation from a control run with no surface 
forcing. The initial front is centered along y = 192 km. 

The numerical simulations presented here treat three processes which were not considered in 
Nofs (1983) original work. Aside from the effects of an along-stream cooling gradient, the 
numerical simulations also consider: 1) the forced meridional response due to a cross-stream 
cooling gradient as outlined in Section III; 2) the effects of momentum mixing during convective 
adjustment; and 3) the time-dependent response of the zonal flow. These three effects are critical 
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in the simulations since they determine much of the response in the immediate vicinity of the 
surface front. 

A. EFFECT OF VERTICAL MIXING OF HEAT ONLY (EXPERIMENT 4-1) 

The surface v components at hour 72 for the simulation with no momentum mixing during 
convective adjustment are shown in Fig. 4-1. Within the region where there is a uniform along- 
stream cooling gradient (see Fig. 2-2), the v components are negative (southward) with maximum 
current speeds near 0.6 cm s’^. The cross-stream components south of the front are in thermal 
wind balance as in in Nof (1983). The magnitude of the v components in this simulation are five 

times smaller than Nofs predictions. Nof used an observed density increase (Gorshkov, 1978) of 

-3-3 . -2 

1 x 10’ gm cm’ along the Gulf Stream extension. A surface cooling rate of 2000 W m’ would 

have to persist for a period of 10 days to obtain such a density increase in the mixed layer. From 

the observations by Agee and Howley (1976) (Fig. 2-4) and Kondo (1977) (Fig. 2-6), such large 

cooling rates appear to be unrealistic, which implies the density increase used by Nof is not due to 

surface cooling alone. 

There is also a relative maximum in the southward velocity immediately to the south of the 
initial zonal current maximum. The process for the inducement of cross-stream velocities near a 
jet in the presence of cross-stream cooling is described in Section III. The magnitude of the 
maximum cross-stream flow in Fig. 4-1 is only slightly larger than that from the two-dimensional 
(y-z) simulations. The slight increase is due to the adjustment of the flow to a weak zonal 
temperature gradient induced by the zonal cooling gradient. Most of the cross-stream velocity in 
the immediate vicinity of the current maximum appears to be due to the two-dimensional effect of 
cross-stream cooling and not the imposed along-stream cooling gradient. At the current 
maximum, the cross-stream flow is about half the magnitude of the cross-stream flow further to 
the south. Thus, displacements of the surface current system would be less than might be inferred 
from the southward velocities that develop where the zonal gradient in surface cooling is a 
maximum. 
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Figure 4-1. Surface v component at hour 72 for Experiment 4-1. Negative values indicate 
southward flow. 



Relatively large zonal gradients in v occur where the zonal gradient of the cooling changes (x 
= 320 km and x = 1000 km), especially in the extreme southern portion of the domain where 
the zonal velocities are small. Between y = 100 and 190 km, the zonal gradient in v is more 
diffuse. Zonal advection of temperature tends to weaken the zonal temperature gradients near x = 
320 km and x = 1000 km. Thus, the southward components, which are in thermal wind balance, 
are smaller in these areas. 

The decreases in the surface temperature field and the change in the surface flow field at 72 
hours are shown in Fig. 4-2. The pattern of the surface temperature changes is very similar to the 
pattern of the surface cooling (Fig. 2-2), which indicates that the changes in the surface 
temperature are due to the surface cooling and not the internal response of the ocean. The 
temperature changes near the western boundary of the surface cooling are slightly modified by the 
tendency of zonal advection to restore the temperature field to the initial conditions. As in the 
case with the velocity components, the maximum temperature change of 0.8°C is very close to the 
value in two-dimensional simulations. 

The changes in the surface flow field are consistent with geostrophic adjustments in response 
to the changes in the surface temperature field. Because the cross-stream cooling gradient reduces 
the surface horizontal temperature gradient, the surface zonal components are diminished as 
denoted by the westward pointing arrows in Fig. 4-2. The maximum reduction in the surface 
zonal components of 8.0 cm s~* is similar to the maximum decrease in the two-dimensional 
simulations and occurs along the current maximum at the downstream boundary of the domain. 
Smaller decreases in the zonal velocity occur to the north and to the south of the front where the 
cross-stream cooling gradient is smaller. West of 600 km, very small decreases in the zonal 
velocity occur due to the combination of a reduced cross-stream cooling gradient and the 
tendency of zonal advection to restore both the temperature and zonal velocity to their initial 
conditions. 
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The contour interval is 0.1 



An 18-day integration using a cooling rate six times smaller than the previous experiment was 
performed to study the effect of having the same total cooling, but over a longer period of time. 
The magnitude of the cooling is very similar to the maximum average cooling during February 
computed by Gorshkov (1978) (Fig. 2-5). The cross-stream velocities at day 18 for the simulation 
with a lower cooling rate are shown in Fig. 4-3. In the immediate vicinity of the front, the cross- 
stream velocities are much smaller than in the previous simulation (Fig. 4-1). This is consistent 
with the analysis in Section III which shows that the cross-stream velocity response due to a 
cross-stream cooling gradient is only weakly time dependent and will be reduced when the cross- 
stream cooling gradient is reduced. The smaller cross-stream components also support the 
conclusion that it is the value of the cross-stream cooling gradient and not the along-stream 
cooling gradient that determines the cross-stream flow in the immediate vicinity of the front. 
South of the front, the meridional flow is slightly smaller than in the previous simulation. In this 
region, the tendency for the along-stream cooling gradient to establish a zonal temperature 
gradient is opposed by the tendency for zonal advection to restore the temperature field to the 
initial condition. The effect is most noticeable where the initial current is strong, between y = 110 
and 190 km. In the extreme southern portion of the domain (y < 110 km) where the zonal 
advective effects are smaller, the meridional flow is comparable in magnitude to the flow in the 
earlier simulation. 

In summary, the combined effects of a cross-stream and an along-stream cooling gradient 
induce a southward flow in these numerical simulations, but the magnitude of that flow is too 
small to produce displacements of the Gulf Stream of the magnitude suggested by Nof (1983). 
South of the current maximum the southward velocities are in thermal wind balance with the 
induced zonal temperature gradient, which is consistent with NoFs (1983) predictions. Near the 
front, the tendency for the formation of a southward flow is due more to the constant adjustment 
of the flow fields to the reduced meridional temperature gradient than the along-stream cooling 
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Figure 4-3. Surface v component at day 18 for a simulation using a cooling rate six 
times smaller than the first simulation of Experiment 4-1. The contour 



gradient. Much of the change in the magnitude in the initial jet can be explained by the reduced 
horizontal temperature (pressure) gradient. A strong short period of cooling is more effective in 
inducing cross-stream velocities in the immediate vicinity the front than is a steady prolonged, but 
smaller, cooling rate. The stronger short cooling event is also more effective in producing larger 
meridional velocities south of the current maximum. Because there is no surface forcing on the 
upstream flow in these simulations, there is a tendency for zonal advection to restore the fields to 
their initial condition. The tendency for zonal advection to restore the fields to their initial 
condition is larger over relatively longer integrations. 

B. EFFECT OF MOMENTUM MIXING (EXPERIMENT 4-2) 

It was shown in Section III that the inclusion of momentum mixing in the convective 
adjustment is critical in determining the magnitude and direction of the response of an intense 
oceanic flow to cross-stream cooling. With momentum mixing, surface cooling causes the 
relatively high zonal momentum near the surface to be mixed with relatively low zonal 
momentum immediately below, so that the surface currents decrease relative to their initial values. 
The horizontal gradient in the surface cooling, and the associated decrease in pressure gradient is 
not large enough to balance geostrophically the decrease in the momentum. Consequently, the 
surface flow becomes sub-geostrophic, and positive (northward) v components occur near the 
front. This behavior is opposite of that predicted (Experiment 4-1) when momentum mixing is 
not allowed. 

The surface meridional velocities at hour 72 for the simulation that includes convective mixing 
of momentum are shown in Fig. 4-4. The maximum cross-stream speeds near the current 
maximum are now near 2.5 cm s as opposed to -0.3 cm s’ ^ when momentum mixing is not 
included. South of the current maximum, the southward flow is only slightly smaller than in 
Experiment 4-1. The flows in this region are similar in the two experiments because there is little 
vertical current shear at the initial time, so momentum mixing has almost no effect. Thus, the 
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the flow is negative (southward). 



drastic effect of including or excluding momentum mixing is confined to the region where the 
initial current has a sizeable vertical shear, i.e. near the front. 

The maximum cross-stream component at the front is similar to the two-dimensional 
simulations, but the contribution from the along-stream differences is also noticeable. At the 
front, the cross-stream velocities near the western boundary of the cooling are only half the 
magnitude of the cross-stream components farther downstream. Between x = 600 km and x = 
1000 km, the surface cross-stream components near the front are nearly constant in x even though 
the forcing is a function of x. The cross-stream response at the front depends more on the initial 
shear in the zonal velocity than the imposed surface cooling in these simulations. The surface 
cross-stream components at the front between x = 300 km and x = 600 km are smaller due to 
the tendency of zonal advection to restore the fields to their initial condition. Also, at the 
southern edge of the front near y = 170 km, there is an area of divergence, which also occurred in 
the two-dimensional simulations. However, the divergence in the three-dimensional simulations is 
larger due to the added effect of the southward meridional flow forced by the along-stream cooling 
gradient south of the front. The divergence, and hence upwelling, is nearly 20% larger in the 
three-dimensional simulation. 

The disturbance with a wavelength of about 100 km along the front in Fig. 4-4 is not due to 
contamination by inertial components that are not removed by the filtering. A time series of the 
instantaneous v component and the filtered v component at x = 1100 km and y = 188 km is 
shown in Fig. 4-5. An inertial oscillation in the instantaneous flow (dashed line) is very evident. 
Although there is a 6 cm s'^ oscillation in the instantaneous v components, the low-pass filtered 
data show no sign of the inertial oscillation. Rather, the 100 km wave is due to the growth of 
perturbations on a baroclinically unstable flow. The wavelength for maximum growth rate for this 
flow was determined from separate simulations (not shown) to be on the order of 75-100 km. 
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Figure 4-5. Time series of the instantaneous v component at the surface (dashed line) 
and filtered v component (solid line) at x = 1100 km and y = 188 km 
for Experiment 4-2. 



The maximum reductions in the zonal components (Fig. 4-6) near the surface jet are almost 
twice as large as in Experiment 4-1. The change in the surface zonal components is due to a 
combination of the adjustment of the flow to a weakened meridional temperature gradient and the 
reduction of the flow through convective mixing with the sub-surface flow. The reductions in the 
zonal components extend farther to the south in this simulation than in Experiment 4-1. The 
larger southward extent of the reductions in the zonal components is due to the divergence of the 
cross-stream flow south of the front, which tends to weaken the horizontal temperature gradient 
there. The effect is most noticeable at the eastern boundary of the cooling where there are 
reductions in the zonal components south of the current maximum, and significantly smaller 
reductions north of the maximum. Away from areas of meridional cooling gradients and areas of 
vertical shear, there is little change in the surface zonal flow. 

The change in surface temperature after 72 hours (Fig. 4-7) is similar to the surface 
temperature changes in Experiment 4-1 with only minor differences. As in Experiment 4-1, the 
temperature changes are primarily due to the surface heat flux rather than the ocean's internal 
response to the forcing. The maximum temperature decreases in Experiment 4-2 are shifted 
slightly to the north and the west of the decreases that occurred in Experiment 4-1. The 
northward shift is due to advection by the positive v components at the front, and the westward 
shift is due to the weakened current maximum and thus, weakened tendency for zonal advection 
to restore the fields to their initial condition. 

The changes in the vector flow field after 72 hours (Fig. 4-7) are also similar to Experiment 
4-1, except in the immediate vicinity of the front. There is a noticeable tendency for the vector 
flow to cross the lines of constant temperature change, indicating that the change in the flow is 
not as geostrophic as in Experiment 4-1. Away from the immediate vicinity of the surface front, 
the changes in the flow are parallel to the changes in the temperature. 
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Figure 4-6. Change in surface u component of velocity at hour 72 for Experiment 4-2. 

Negative values indicate a reduction in the eastward current maximum. 

The contour interval is 2 cm s“l . 
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In summary, the cross-stream components in the immediate vicinity of an ocean front are 
larger, and in the opposite direction, to the simulation with no momentum mixing during 
convective adjustment. South of the current maximum, in the region of weak vertical shear, the 
induced southward flow is only slightly changed by the inclusion of vertical mixing of 
momentum. The reductions in the zonal component near the current maximum are twice as large 
as those that occurred in Experiment 4-1. As in Experiment 4-1, the changes in temperature are 
primarily due to the surface forcing rather than an internal adjustment of the ocean. The effects of 
vertical mixing of momentum dominate the effects of the imposed surface cooling gradients near 
the current maximum. The tendency to produce northward v components at the front when 
momentum mixing is included is large enough to compensate a cross-stream flow of the 
magnitude predicted by Nof (1983). The above conclusions obvious^ depend on the assumed 
initial vertical profile of zonal momentum, and on the assumed form of the vertical mixing of 
momentum (adjustment based on the value of the local Richardson number). At least for these 
simulations, strong cooling events alone cannot explain the observed southward shift of the mean 
position of the Gulf Stream during winter. 

C. EFFECT OF WIND FORCING (EXPERIMENT 4-3) 

It was suggested in Section III that an increase in the zonal wind stress during winter may be 
an alternate explanation for the observed southward shift of the Gulf Stream. An additional 
integration was performed to consider the response of a strong ocean current as it flows into a 
region of wind forcing. The wind stress increases in the along stream direction from 0 at x = 320 
km to 0.2 N-m at x = 1000 km. There is no meridional dependence of the wind stress. To 
isolate first the effects due to wind forcing, there is no surface cooling specified in Experiment 4-3. 
The choice of wind stress which increases in the along-stream direction is supported by 
observations by Agee and Howley (1976) (Fig. 2-4). Simulations with a meridionally dependent 
wind stress were performed but did not significantly alter the results presented in this section. 
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The changes in the surface temperature field after 72 hours are shown in Fig. 4-8. The 
reductions along the northern boundary are due to the effect of upwelling forced by an Ekman 
divergence and are not of interest for this study. Reductions in the sea-surface temperature are 
evident in the eastern two-thirds of the domain, which suggests that the front has been shifted 
southward consistent with the direction of the Ekman flow. The maximum decreases along the 
current maximum are in excess of 1.4°C at the eastern boundary of the domain, which is 
equivalent to a 9 km southward shift in the front. The shifts in the previous simulations, which 
are due to cooling alone, are less than 1 km. The 9 km shift is comparable to the shift noted in 
the two-dimensional simulations. The temperature decreases farther to the west are not as large 
as those in the eastern domain due to the smaller southward advection where the wind stress is 
smaller. This integration was extended to 7 days to determine a limit for the southward shift of 
the front for this wind forcing. After seven days of forcing, the front was displaced 13 km 
southward. The rate of the southward advection of the front has been slowed which indicates a 
limit on the effectiveness of wind stress for shifting the front. As the front is advected southward, 
the slope of the front becomes larger, and the depth to which the frontal region is well-mixed 
becomes greater. Consequently, the effects due to advection become smaller since the Ekman 
(advecting) velocities are inversely proportional to the mixed-layer depth. 

An integration similar to Experiment 4-3 but also including cooling similar to that described in 
Experiments 4-1 and 4-2 was performed to investigate the combined effects of the atmospheric 
forcings on strong ocean currents. The changes in sea-surface temperatures after 72 hours (not 
shown) are similar to the pattern in the wind-forcing only experiment, but the decreases in surface 
temperature are larger (especially south of the front) due to the extra effect of surface cooling. The 
effect of the cross-stream components (not shown) is also additive. The meridional flow south of 
the current maximum is enhanced due to the adjustment of the flow to the forced along-stream 
temperature gradient. At the front, the meridional flow is still southward, but at a reduced speed 
relative to the case with no surface cooling. 
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As in the two-dimensional simulations, the effect of wind forcing is more effective in shifting 
the position of a intense current system than are very strong gradients in the horizontal cooling. 
However, none of the simulations support the hypothesis that atmospheric forcing is capable of 
shifting an intense oceanic baroclinic zone (and the associated current system) far enough to cause 
significant feedback to the atmosphere. To this point, the effects of atmospheric forcing on a 
uniform zonal flow only has been considered. The Gulf Stream is a meandering current system 
that adjusts to changes in the bottom topography and perturbations on a strong baroclinic flow. 
It is possible that atmospheric forcing may affect the processes that produce meanders in the Gulf 
Stream, and ultimately lead to changes in the mean path of the flow. 
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V. THE EFFECTS DUE TO TOPOGRAPHY 



The position of the Gulf Stream has been observed to change rapidly as it encounters 
topographic features associated with the New England Seamount Chain (NESC) (Vastano and 
Warren, 1975). Richardson (1981) observed that the NESC appears to be responsible for a 100 
km southeastward deflection of the Gulf Stream and for a quasi- permanent warm ring-meander to 
the south of the seamount chain. Richardson also observed that areas downstream of the NESC 
are favored locations for large-scale meander development. The seamount chain affects the Gulf 
Stream because of the height of its topographic features, which in places rise to 3 km above the 
ocean floor, and also because the flow associated with the Gulf Stream penetrates deep into the 
water column. 

In this section, the effect of atmospheric forcing on baroclinic flow moving over topography is 
studied to determine if the flow can be steered away or toward a topographic feature. A second 
objective of this section is to determine if the atmospheric forcing affects the actual adjustment of 
the flow to underlying topography. If the adjustment of the flow is sensitive to the surface forcing, 
then large changes in the position of the current system may occur. If the flow can be steered 
away from or toward topographic features, then the character of the flow downstream of the 
topography may change radically as observed by Richardson (1981). In either case, it may be 
possible to affect significantly the position of the surface front over time scales that are important 
to the atmosphere. 

The simulations described in this section are very similar to a linear study by Smith (1984) 
who looked at the two-dimensional response of atmospheric flow over a mountain. In that study, 
as well as in this, the direction of the flow reverses with height. The reversal of the current with 
height in this study is due to the assumption of no barotropic component. The current reversal in 
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Smith's study was a necessary condition for the propagation of disturbances away from a 
topographic feature. Although the purpose of this section is not to reproduce Smith's results, his 
analytic work does provide a guideline for determining if the flow is adjusting to the topography 
in a realistic fashion. 

An example of a cross-section of the initial zonal velocities used in this study is shown in Fig. 
5-1. The flow is eastward in the upper 400-500 m and westward in the lower layers. Thus, the 
flow in the layer adjacent to the bottom is from east to west. The topography used in these 
simulations is a Gaussian-shaped seamount that rises 200 m above a 1 km deep ocean floor. The 
width of the seamount is approximately 120 km. Although this investigation considers the effect 
on stratified flow only, the conservation of potential vorticity is expected to be very important in 
these simulations as in Smith's (1984) results. The initial conditions are similar to the those in 
Section IV. The surface front is centered along y = 128 km for all the simulations, and the initial 
zonal velocities are prescribed to be in geostrophic balance 

u = 4 (P v ~ 7T p * • t 5 - 1 ) 

P o 1 D 

These simulations of the adjustment of the flow to topography neglect the barotropic 
component. In an atmospheric study of lee cyclogenesis, Hayes (1985) incorporated topographic 
effects into Eady's (1949) solution for baroclinic growth. Hayes found that the steady state 
solution for flow over a mountain decayed exponentially with height above the mountain. The 
vertical structure in these ocean simulations can be estimated (Fig. 5-2) from Hayes' solution 
using values of 100 km for the mountain width, 10 km for the baroclinic radius of deformation 
and -0.5 for the ratio of the mean current at the botttom to the change in the mean current up to 
the thermocline. Notice that the amplitude of the response changes sign with height and is almost 
totally baroclinic. The barotropic component does not seem to be important for the adjustment 
of the flow for the topographic scales in this study. Hayes shows that the vertical extent of the 
topographic effect is reduced as the horizontal scale of the topography decreases relative to the 
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Figure 5-1. South-to-north cross-section of the initial zonal velocities used in the 
simulations which include effects due to topography. The contour interval 
is 20 cm s _1 , and positive flow is toward the east. 
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Figure 5-2. Dependence of the amplitude of the response of the flow on non-dimensional 
height for the initial conditions and topography used throughout 
Section V • 



deformation radius. Since the wavelength of the seamount in the present study is relatively long 
compared to the baroclinic radius of deformation, the topographic effect should be felt all the up 
to the surface. However, as shown in Fig. 5-2, the response will be primarily baroclinic as 
assumed in the numerical experiments. 

The investigation is divided into two sets of experiments. The first experiment considers the 
effect on the flow when the maximum surface current will be forced away (south) from the 
seamount, and the second experiment considers the response of the flow when the maximum 
surface current is being forced toward the highest point of the seamount. Each experiment 
contains three simulations. The first simulation of each experiment is a control case which is 
integrated to 6 days to determine the response of the flow to topography alone. The second 
simulation is identical to the first except that the topography is shifted 16 km north relative to the 
first simulation. Thus, the first two simulations of each experiment highlight the adjustments of 
identical flows which are initially over a slightly different topography. The third simulation uses 
the fields from day 3 of the first simulation and then includes an eastward wind stress which is 
applied impulsively. The tendency of the Ekman advection in this simulation will be to shift the 
front south as in Sections III and IV. Surface wind stress is chosen to force the problem since it 
was shown to be effective in shifting the position of the surface current system in Sections III and 
IV. The wind stress increases from zero at x = 120 km to a maximum value of 0.2 N-m at x = 
260 km. The wind stress has no meridional dependence. 

By comparing the third simulation with the first simulation, it is possible to determine how the 
surface forcing affects the adjustment of the flow to the topography. By comparing the third 
simulation with the second simulation, it is possible to determine if the atmospheric forcing is 
effective in steering the flow toward or away from the topography. If the surface forcing is able to 
steer the flow, the adjusted flow in third simulation will resemble the second simulation more that 
the first. If the surface forcing is not effective in steering the flow, the adjusted flow will resemble 
the first simulation. 
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A. SURFACE CURRENT FORCED AWAY FROM SEAMOUNT (EXPERIMENT 5-1) 

The first simulation considers the response of flow over a seamount which is centered at x = 
640 km and y = 128 km. In this case, the maximum surface flow is directly over the highest 
point of the seamount. The surface temperature at day 6 for the simulation with no surface 
forcing is shown in Fig. 5-3. Upstream of the seamount near x = 600 km, the fluid particles 
encounter a shallower water column and curve anticyclonically consistent with the conservation of 
potential vorticity. Downstream of the seamount, the isotherms curve cyclonically as the flow 
moves over deeper water. Without surface forcing, the front is approximately at the same latitude 
downstream of the seamount as it is upstream of the seamount. 

The conservation of potential vorticity is a powerful dynamic constraint in this simulation. 
However, because the fluid is stratified, a discussion in terms of density can explain the curvature 
in the surface isotherms. Fig. 5-4 is a west-to-east cross-section of the temperature field at day 6. 
A dramatic change from the initial temperature field is evident (there is no zonal dependence in 
the initial conditions; the isotherms are initially flat). There is a marked lifting of the isotherms on 
the eastern side of the seamount, while on the western side, the isotherms are depressed. Since the 
shear flow near the bottom is from east to west (Fig. 5-1), there is convergence (divergence) on 
the eastern (western) side of the seamount which induces an upward (downward) vertical motion. 
The convergence (divergence) at depth is balanced by a divergence (convergence) above so that 
the net divergence over the water column is zero. Notice that the upwelling on the eastern side of 
the seamount penetrates the thermocline. The pattern of vertical motion produces a column of 
relatively cold water above the eastern side of the seamount and relatively warm water above the 
western side. The zonal pressure gradient induced by this warm-cold pattern is negative above the 
western side of the seamount, positive over the seamount and negative above the eastern flank of 
the seamount. The surface meridional geostrophic flow for these pressure gradients is northward 
above the eastern and western sides of the seamount and southward over the seamount. Thus, v^ 
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Figure 5-3. Surface temperature at day 6 for the first simulation of Experiment 5-1. 
The contour interval is 1.0 °C. 
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Figure 5-4. West-to-east cross-section of temperature across y = 128 km at day 6 for 
the first simulation of Experiment 5-1. The contour interval is 1.0 °C. 



< 0 (anticyclonic) on the western slope of the seamount and v x > 0 (cyclonic) on the eastern 
slope of the seamount. In this simulation, potential vorticity is conserved by an adjustment to 
the local areas of upwelling and downwelling generated as the flow moves over the topography. 
The train of lee waves to the west of the seamount is an effect that was noticed by Smith (1984). 
The appearance of these waves in this study is an indication that the flow is adjusting in a fashion 
similar to previous investigations. 

The second simulation of this experiment considers the response of flow over a seamount 
which is centered at x = 640 km and y = 144 km. Thus, the initial surface flow in this 
simulation is south of the center of the seamount. The surface isotherms at day 6 are shown in 
Fig. 5-5. The response of the flow is very different from the first simulation. Although the initial 
flow is located over smaller topographic features, the amplitude of the meander is nearly two 
times larger than in the earlier simulation. The larger amplitude meander develops because the 
initial adjustment of the surface flow west of the seamount is anticyclonic to the north and the 
flow is forced toward larger changes in the topography. By contrast, the flow in the first 
simulation adjusts so that it is shifted away from the largest slopes in the topography. 

The first simulation is repeated with a 0.2 N-m wind stress in the eastward direction which is 
applied impulsively at day 3. The difference in the surface temperatures at day 6 between the 
simulations with and without forcing is shown in Fig. 5-6. The decreases in the surface 
temperature are due to a southward advection of the front. The maximum decrease of 2.0 °C 
occurs along the current maximum and is comparable to the decrease described in Section IV for 
a flat bottom case. The contours of temperature difference parallel the surface flow which 
indicates a uniform advection. Along the northern boundary of the domain, the decreases are due 
to upwelling forced by an Ekman divergence at the boundary and are not of interest in this study. 

The temperature response below the surface layer (not shown) is remarkably similar to that in 
the simulation without surface forcing. The fluid below the surface layer is protected against the 
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Figure 5-5. As in Fig. 5-3 except that the seamount is centered at x = 640 km and 
y = 144 km. 
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Figure 5-6. Difference in surface temperature fields at day 6 for Experiment 5-1 

between the simulation which is forced by an eastward wind stress and the 
simulation which is not forced. The contour interval is 0.5 °C. 



action of the wind stress by the initial stratification. Thus, the flow adjusts in a similar manner as 
when forcing is not included, and does not depend on the upper layers in this simulation. When 
the simulations are extended to 8 days, the effect of Ekman advection becomes reduced as the 
front steepens, as in Experiment 4-3. The surface temperature decreases are only slightly larger 
(2.5 °C) than at day 6 and still parallel the surface flow. 

Although Ekman advection has shifted the surface flow south of the seamount in the wind- 
forced case, the flow does not develop as large a meander as in the second simulation when the 
flow is initially south of the seamount. Thus, the forcing does not affect the actual adjustment of 
the flow in these simulations. Also, the forcing has not steered the flow so that the adjustment 
evolves from the smaller amplitude meander of the first simulation to the larger amplitude 
meander of the second simulation. The meander in the third simulation is similar to the meander 
in the first simulation. 

B. SURFACE CURRENT FORCED TOWARD SEAMOUNT (EXPERIMENT 5-2) 

The three simulations of Experiment 5- 1 were repeated with the highest point of the seamount 
shifted southward 40 km for each case. In the these simulations, the initial surface current 
maximum is always north of the center of the seamount, which is along y = 88 km in the first 
simulation. The surface isotherms for the control simulation without forcing at day 6 are shown 
in Fig. 5-7. The amplitude of the meander is much larger in this simulation (100 km) than in any 
of the simulations of Experiment 5-1. The amplitude of the meander is large in this simulation 
because the adjustment of the flow to the topography forces the flow over the highest point of the 
seamount. Although the amplitude of the meander is about 100 km, the surface front is 
approximately at the same latitude west of the seamount as it is east of the seamount. 

A 16 km shift in the initial conditions produces a radically different response in the amplitude 
of the meander. The simulation with the seamount centered at x = 640 km and y = 104 km 
produces a surface isotherm pattern that is very similar to the pattern in the first simulation of 
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Figure 5-7. As in Fig. 5-3 except for Experiment 5- 



Experiment 5-1. (Fig. 5*3). As the adjustment of the flow to the topography is very sensitive at 
this location, it may also be sensitive to the effects of surface forcing. 

The first simulation of this experiment is repeated with the effect of wind forcing included to 
determine whether or not the forcing is able to shift the adjustment of the flow from the large 
amplitude meander to the smaller amplitude meander. The differences in the surface isotherms at 
day 6 between the integration which included wind forcing and the first simulation are shown in 
Fig. 5-8. Notice that the temperature differences in this experiment tend to parallel the surface 
flow as in Experiment 5-1. Also, the amplitude of the meander has not changed. Although Ekman 
advection has shifted the surface front 10 km in this simulation, the effect of the wind stress is 
confined to the upper layers by the initial stratification, and the adjustment of the flow does not 
change from the large amplitude meander mode to the smaller amplitude meander mode. 

In summary, the response of the flow in the simulations without surface forcing is consistent 
with the conservation of potential vorticity. When the effects of an eastward wind stress are 
included, Ekman advection induces a uniform southward displacement of the surface front. 
Although the position of the surface current system is changed due to advection, the adjustment 
of the sub-surface flow to the topography and the amplitude of the meanders which are produced 
does not change in these simulations. The effects due to atmospheric forcing cause no significant 
changes in the position of the front or changes in the character of the downstream flow in these 
simulations. 
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Figure 5-8. As in Fig. 5-6 except for Experiment 5-2. 



VI. SUMMARY AND CONCLUSIONS 



The objective of this study is to investigate the effect of atmospheric forcing on the location of 
the mean path of the Gulf Stream. Associated with the Gulf Stream is a strong baroclinic zone 
that extends eastward from the continental United States for hundreds of kilometers. A shift in 
the position of the Gulf Stream may have important consequences for the atmosphere such as in 
the case for predicting explosive cyclogenesis events (Sanders and Gyakum, 1980). One of the 
emphases of this study has been on the role of surface cooling in determining the path of the Gulf 
Stream. Nof (1983) has proposed that an observed southward shift of the Gulf Stream of 90 km 
during winter may be due to the effects of surface cooling. This investigation deals mainly with 
the ocean's response to strong atmospheric forcing on time scales of about 3 days, which should 
contribute to the seasonal effects of surface forcing. 

The two-dimensional response of the current system to surface cooling is highly dependent on 
the vertical mixing of both heat and momentum in these numerical simulations. In a convective 
overturning case with no momentum mixing, the magnitude of the horizontal pressure gradient is 
reduced, which induces a weak cross-stream flow toward the warmer water. A thermally indirect 
circulation is centered on the warm edge of the surface front. Although the sense of the cross- 
stream current is to advect the front southward at the surface, the vertical circulation acts to 
increase the horizontal temperature gradient and maintain the front at its initial position. When 
momentum mixing is also included, the decrease in the along-stream current at the surface is 
larger than the change in the horizontal pressure gradient, and a cross-stream velocity toward the 
cooler water occurs near the surface. The v components are an order of magnitude larger than in 
the simulation that did not include the momentum mixing, but they are in the opposite direction. 
The sense of the vertical circulation is also changed when the momentum mixing is allowed. The 
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vertical circulation becomes thermally direct and tends to diminish the horizontal temperature 
gradient. As high momentum is transferred downward into a layer of diminished horizontal 
pressure gradient, strong cores of southward velocity occur immediately below the surface layers. 

During winter, there is an increase in the mean eastward wind stress associated with the 
stronger and more frequent synoptic disturbances in the atmosphere. The inclusion of a 
southward wind stress to simulate a cold outbreak from the north in the two-dimensional model 
produces an Ekman-type advection above the main thermocline. The Ekman component is 
directed upstream which reduces the magnitude of the initial jet. The magnitude of the cross- 
stream velocity is also reduced relative to the simulation with no wind stress. Except for a slight 
southward shift, the net effect on the temperature field is minimal. Imposing an eastward wind 
stress in the two-dimensional model to simulate a cold outbreak from the west, induces a larger 
response than does the southward wind stress. The sense of the Ekman advection is to displace 
the front southward which affects both the along-front velocities and the temperature field. Due to 
the strength of the Ekman advection, the cross-stream velocities are in the opposite direction from 
the simulation without the wind stress. The divergence at the front and the horizontal gradient in 
cooling both act to reduce the intensity of the surface horizontal temperature gradient. 

In these simulations, a modest wind stress of 0.2 N-m has a greater effect on the position of 
the Gulf Stream than does a very strong horizontal gradient in the cooling. Ekman advection due 
to an eastward wind stress displaces the current system southward and increases the vertical tilt of 
the front. The strong response in the simulation with an eastward wind stress suggests an alternate 
explanation for the observed shift in the Gulf Stream position and increased vertical tilt of the 
associated temperature front (Worthington, 1976). Horton (1984) also found that the main 
response of the Gulf Stream to forcing by hurricane Dennis was that the surface layer was shifted 
by Ekman advection. 
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The three-dimensional response of an intense oceanic current system as it enters a region of 
cooling was also numerically simulated. In the immediate vicinity of the associated surface front, 
the flow response is primarily due to a cross-stream gradient in cooling rather than the along- 
stream cooling gradient. As in the two-dimensional simulations, relatively small southward 
velocities (-0.3 cm s'*) occur when momentum is not mixed vertically during convective 
adjustment. When momentum mixing is included, larger northward velocities (2.5 cm s’*) occur. 
Performing an 18-day integration with a reduced but steady cooling function produced a very 
weak (-0.1 cm s’*) cross-stream response in the immediate vicinity of the front. These simulations 
do not support the hypothesis that downstream surface cooling alone is responsible for an 
observed southward shift in the mean position of the Gulf Stream during winter. 

Away from the immediate vicinity of the front, the response.in the flow field is entirely due to 
the prescribed along-stream cooling gradient. In both experiments, a southward flow field 
develops in response to cooling that increases toward the east. The cooling produces a zonal 
pressure gradient to which the flow adjusts geostrophically. That is, the meridional component is 
in the thermal wind balance as in Nof (1983). 

Because the upstream flow is not forced, there is a strong tendency in the three-dimensional 
simulations for zonal advection to return the temperature and flow fields to their initial condition. 
The effect is most evident near the region where the strong currents first encounter cooling. The 
restoring advection is not important in areas of weak zonal flow, and effects due to the cooling 
can be seen farther upstream than in areas where the initial zonal flow is strong. 

The inclusion of a moderate increase in the eastward wind stress shifts the surface front 
southward 9 km after 3 days of integration. However, an additional 4 days of integration results in 
a further shift in the surface front of only 4 km. The steepening of the frontal interface limits the 
effectiveness of the wind stress in advecting the front farther southward. 
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Investigations into the large-scale forcing of these current systems to explain observed seasonal 
shifts seems warranted. The position of the Gulf Stream is not displaced enough by the local 
atmospheric forcing to significantly affect feedback by the ocean to the atmosphere on the time 
scale of 3-6 days. Although an increase in the zonal wind stress is able to displace the upper layer 
of a strong current system in these simulations, the effect does not extend below the thermo cline. 
The sub-surface front is a effective anchor of the current system in these simulations. 

Because the Gulf Stream penetrates very deep into the water column, the effect of gradients in 
the underlying topography may be important in determining the path of the Gulf Stream. A 
primitive equation model in a-coordinates is used to simulate the adjustment of ocean flow to 
bottom topography. The conservation of potential vorticity is a strong dynamical constraint for 
flow that is adjusting to topography in the numerical simulations. The flow acquires anticyclonic 
curvature as it moves toward the seamount, and cyclonic curvature as it move away from the 
seamount. The inclusion of a surface wind stress produces a nearly uniform Ekman advection of 
the surface front, but does not affect the flow's adjustment to the underlying topography. The 
initial stratification protects the lower layers against the effects of the wind and from affecting the 
main adjustment of the flow. 

The experiment was repeated with the seamount moved south 40 km relative to the maximum 
surface current. In control experiments, the maximum surface current has a much larger 
amplitude meander than occurred when the initial current maximum is directly over the 
seamount. The adjustment of the flow to the topography shifts the front so that the maximum 
currents flow over the highest points in the topography and larger amplitude meanders occur. The 
inclusion of surface forcing is not able to shift the adjustment from the large amplitude meander 
to the smaller meander. The difference in the surface temperatures is parallel to the surface flow, 
though not as uniformly distributed along the flow as in Experiment 5-1. The inclusion of 
atmospheric forcing does not affect the adjustment of the flow to the topography so that there 
will be feedback to the atmosphere on time scales of 3-5 days. 
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This study numerically simulates the effects of atmospheric forcing on intense ocean currents. 
These simulations provide little support for the hypothesis that surface forcing is capable of 
shifting ocean current systems so that there is an effect felt by the atmosphere over time scales of 
about three days. However, each of the experiments isolates effects which are important to the 
local balance of forces near strong ocean currents. These simulations aid in understanding the 
response of the ocean on short-time scales, and hopefully will be helpful in the understanding and 
prediction of the large-scale flow. 
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APPENDIX A 



STABILITY ANALYSIS 

The instantaneous flow fields along the current maximum display an oscillatory behavior in 
the experiments when momentum mixing is included. These oscillations have a period near 15 
hours rather than the 21 -hour inertial period. The purpose of this appendix is to demonstrate the 
effect of advection on the inertial frequency in a numerical model by use of a stability analysis. 

Consider the equations for a flow where changes in the velocity are due to inertial motion and 
zonal advection only. The linearized equation for the horizontal flow is 

w t = — U w x — ifw (A.l) 

where w = u + iv. Numerical simulations involve approximations to the derivatives and may 
alter the frequency of the response. Consider the numerical approximation to (A.l) 

wf +1 _ wj*- 1 = -JM (wf +1 - wf 4 ) - 2i f At wj 1 (A. 2) 

where the subscripts refer to the position in space and the superscripts refer to levels in time. If a 
solution for w is assumed to have the wave form 

wj 1 = W" exp(i k j Ax) 



then one obtains 

r +1 = W"' 1 - 2iW"(i^ sin(kAx) + fAt) (A.3) 

Ax 

Assume an amplification matrix, n, such that W n+ ^ = £i\V n and p Kv n = W n Substitution 
into (A.3) yields a quadratic for p 

f + 2i (-^1 sin(kAx) + f At) n -1 = 0, (A.4) 

Ax 



or solving for /i 
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It = -i (S sin(kAx) + f At) ± (S 2 sin 2 (kAx) - 2 f At S sin(kAx) - f 2 At 2 + l) l/2 (A.5) 

where S = U At/ Ax. 

The radical is always real for the values of S used in the present study. It is possible to obtain 
phase information from the expression for n . The amount of phase change <p in a single time step 
in the numerical solution is given by 

■P = tan-* (A.6) 

real (n) 

For <p to be zero, f At must be equal to -S sin (k Ax). Unlike the analytic solution, <p can never be 
zero for the values of S in this study, because the argument of sin"* would have to be greater than 
one. The period of the response for a given wavelength is given by 

T = — it (A. 7) 

<P 

The oscillations along the current maximum in the numerical simulations have a wavelength 
near 75 km. Using a basic current speed of 2.0 m s‘* and wavelength of 75 km in (A.5), (A.6) and 
(A. 7) yields a period near 15 hours for the oscillations. This result has several implications for the 
simulations presented in this investigation: 1) the equation (A.l) is describing much of the time 
dependence of the flow in the complete numerical model; and 2) zonal advection appears to be 
responsible for the shift in the response away from the inertial frequency. 

Another interesting result of (A.6) and (A. 7) is that the change in the advective-inertial 
frequency is dependent upon the direction of the advecting flow. For U > 0, the period of the 
response will decrease, and for U < 0, the period of the response will increase. The range of the 
periods for different wavelengths at 36°N is between 15 and 21 hours when U > 0, and 21 and 27 
hours when U < 0. 
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